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T— I environment and find an approximately linear dependence on the external field, in accordance 

with lattice data. 
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It is believed that experiments with ultrarelativistic heavy ion collisions can provide access 
to some features of the vacuum of QCD. In this regime the QCD coupling constant is very strong, 
making any investigation, experimentally and theoretically, very involved. Perturbation theory does 
not apply in this case and non-perturbative methods indicate the existence of non-trivial topological 
solutions, accessible only via semi-classical treatments. It was proposed in Ref. [1] that if such 
solutions are indeed present in the strongly interacting matter generated after a heavy ion colli- 
sion, its association with the strong magnetic field produced in non-central collisions can induce a 
phenomenon that has been called "chiral magnetic effect" (CME). This effect generates a current 
of the deconfined quarks that causes a separation of electric charges in the final hadronized state. 
Some suitable observables can be defined in order to analyze the charge separation caused by the 
CME [2]. 

In this work we investigate the effects of a magnetic background on the QCD phase transitions 
[3, 4] using the linear sigma model (LSM) coupled to quarks and to the Polyakov loop (PLSMq). 
The chiral and deconfinement transitions are temperature driven and a finite-temperature mean- 
field approach is implemented in order to model this scenario. 

The confining properties are described by the Polyakov loop L, whose expectation value works 
as an order parameter that indicates whether the system resides in the confined or deconfined phase: 
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where A4 is the matrix- valued temporal component of the Euclidean gauge field and the symbol 
^ denotes path ordering. The integration takes place over the compactified imaginary time T. 

On the other hand, the chiral sector of the theory is described by the LSM part of the model, 
whose degrees of freedom are the a and n mesons. The expectation value of the a meson field 
works as an order pai^ameter for the chiral symmetry, indicating if this symmetry is broken or 
restored, depending on the temperature: 

^, ■ , f (a) 7^ , low r = (a, n) , 

Chiral symmetry : < , Z ) ' ^ l . (2) 

^ ^ [ (a) = , highr ' TT = (7r+,7r°,7r-). 

The spontaneous and explicit chiral symmetry breaking features are encoded in the LSM potential: 

V^{(7,7i) = ^{(j^ + 7i^-v^)^-ha. (3) 

Hot quarks provide a background in which the condensates evolve. The interaction of the quarks 
with the mesons is implemented in a Yukawa fashion, with the following Lagrangian: 

■^q = Wf [ifd^ - g{<^ + '75? • n)] Yf , (4) 

where g is the coupling constant. 

The static uniform magnetic field, pointing in the z direction, is introduced by a replacement 
of the simple derivative by a covariant derivative in the quark kinetic term. This derivative contains 
an Abelian gauge field which encodes all the relevant information about the magnetic field. In the 
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same way, the interaction with the Polyakov loop is introduced via a non-Abelian gauge field in the 
covariant derivative as follows: 



ID = f{d^-iQa^-iAy,), 
and the fuli Lagrangian can now be written as: 



(5) 



^ = W[iS> - gio + irsr • 7?)] 2 [(di^of + (d^K^'f] + jD^Tr+p - ^^(ct, tt) - Vl{L, T) , (6) 

where = iea^ is the covariant derivative acting on coioriess charged pions, represent- 
ing the 4-potentiai for the external magnetic field, {a, n) is the LSM potential and Vl{L, T) the 
Polyakov loop potential, given by [5] : 
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- 270MeV being the critical temperature in the pure SU (3) gauge theory and qq = 
16;r^/45 « 3.51, a\ = —2.47, a2 = 15.2, and bi, = —1.75 are the phenomenological parameters. 
The one-loop corrections to the free energy D. coming from quarks can be written as: 
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The integration measure in the presence of the magnetic field S for T = and T > becomes, 
respectively, as follows: 
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The presence of the magnetic field generates a sum over the Landau levels n in addition to the 
standard finite-temperature sum over the Matsubara frequencies 



(Oi = 2KT£, 

labeled by the integer number £. 

The expectation values can be obtained minimizing the free energy, 

a{o,L;T,B) = V^{o,K) + VL{L,T) + aq{o,L,T), 



(11) 



(12) 



at fixed values of temperature and magnetic field. The interaction piece Q.q{a,L,T) can be sep- 
arated in two parts: a vacuum part - which is temperature independent but still magnetic -field 
dependent - and a thermal piece. Up to one loop the vacuum part is given by: 

det(/ZZ)(''') - m„ 



K 



f=u,d 




I kf\B 



(13) 
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minus the standard vacuum correction in the absence of the magnetic field, 



f=u,d 



'm. 



(14) 



where we have defined the integral 



(15) 



Here capital M means an effective mass which takes into account the effect of the magnetic field B: 



Mlf = m) + ln\qfB\. 



(16) 



The thermal piece is given by: 
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where the dispersion relation for quarks in the external magnetic field is 

(0,n{Pz,o) = [ml{o)+pl + {2n+l-2s)\q\BY^^ , 



(18) 



with constituent quark mass mq{o) = go. The thermal part (17) includes interaction of the Polyakov 
loop with quarks and therefore it breaks the Z3 symmetry, which was initially present in the poten- 
tial for the Polyakov loop Vi. 

The vacuum correction has remarkable consequences on the phase structure. In previous work 
involving the LSM, this connection was usually not taken into account. However, in a theory with 
spontaneous symmetry breaking, the presence of a condensate will modify the masses, so that con- 
tributions from vacuum diagrams can not be subtracted away as trivial zero-point energies. These 
contributions were shown to play an important role at finite temperature by the authors of Ref. 
[6], who incorporate scale effects phenomenologically. Vacuum contributions were considered at 
finite density in the perturbative massive Yukawa model with analytic exact results up to two loops 
in Refs. [7, 8] and, more specifically, in optimized perturbation theory at finite temperature and 
chemical potential in Ref. [9]. More recently, this issue was discussed in different scenarios for the 
chiral and deconfining transitions (see, e.g., Refs. [10]). The exclusion of the vacuum coiTcction 
was motivated by the phenomenological argument that quark degrees of freedom are absent in the 
vacuum. On the other hand, models like the Nambu-Jona-Lasinio (NJL) include such vacuum cor- 
rections automatically. Which treatment represents a better description of QCD is still subject of 
debate, so that we opted for studying both cases. 

For finite temperature and 5 = we obtain a crossover for both transitions, as shown in Fig.l. 
We define the critical temperature as the one where the curve changes curvature, which happens 
simultaneously for the chiral and deconfinement transitions. 
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Figure 1: Expectation values of the order parameters for the chiral and deconfinement transitions as func- 
tions of the temperature. The filled circles represent the CT-condensate and the empty circles stand for the 
expectation value of the Polyakov loop. Both lines are smooth functions of the temperature indicating the 
presence of a crossover. In this plot the condensates are dimensionless. 

At zero temperature, the presence of the magnetic field enhances the breaking of chiral sym- 
metry, making the potential of the chiral field a deeper and increasing the value of the chiral 
condensate (in the zero-temperature case the Polyakov loop variable does not play a role). This 
enhancement effect is known as magnetic catalysis [11], and is illustrated by the left plot of Fig. 
2, where the behavior of the zero-temperature potential for several values of the magnetic field is 
shown. 




Figure 2: Left: effective potential for the condensate a = ^ v at zero temperature for several values of the 
magnetic field B. Magnetic catalysis is manifested in the increase of the condensate and in the deepening of 
the potential. Right: expectation value of the condensate as a function of the magnetic field. The black dots 
are obtained from the PLSMq and the orange Une is a linear fit (19). 

It is interesting to study the dependence of the chiral condensate on the magnitude of the 
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magnetic field. Analytic calculations in chiral perturbation theory [12] and results from lattice 
QCD [13] indicate a linear dependence of the condensate on the strength of the magnetic field. Our 
model is consistent with this picture, and a linear curve fits well the results obtained within our 
approach, as can be seen in the right plot of Fig. 3. The best linear fit in the plot is 

(^ = - = 1.017 + 0.0181^51, v = 87.73 MeV. (19) 

V 

The small deviation of the curve from the linear behavior is qualitatively similar to the one found in 
the context of the Sakai-Sugimoto holographic model [14], although smaller. For a small strength 
of the magnetic field, the dependence of the constituent quark mass (which plays the same role 
as the chiral condensate in our approach) on the magnetic field in Ref. [14] shows a quadratic 
behavior, which can be observed in our plot as well. 

For the full potential, including thermal and magnetic effects, we studied separately the cases 
with and without vacuum coiTcctions. We obtained phase diagrams for both cases, showing that the 
effect of vacuum corrections is indeed remarkable. Fig. 3 shows the corresponding phase diagrams. 



no vacuum corrections with vacuum corrections 




Figure 3: Phase diagram in the B-T plane. Left: without B-dependent vacuum corrections. The critical 
temperatures of the deconfinement (dash-dotted line) and chiral (dashed line) transitions coincide all the 
way, and decrease with B. Right: with Z?-dependent vacuum corrections: the critical temperatures of the 
deconfinement (dash-dotted line) and chiral (dashed line) transitions coincide at B = and split at higher 
values of the magnetic field. A deconfined phase with broken chiral symmetry appears. The vertical line 
represents a typical magnitude of the magnetic field expected to be realized in heavy-ion collisions at the 
LHC [18]. 

If the vacuum corrections are not included, the magnetic field turns both chiral and deconfine- 
ment transitions into first order transitions. Both transitions happen at the same critical temperature, 
which decreases as the magnetic field increases. If included, the vacuum corrections dominate the 
potential, changing drastically the picture. The transitions become crossovers and split as the mag- 
netic field is increased. In both scenarios the presence of the magnetic field have important conse- 
quences, either turning the transitions into first order or causing a split between them, generating a 
new phase. It also reinforces the breaking of the Z3 symmetry that occurs when the Polyakov loop 
interacts with quarks. 

The scenario with vacuum corrections included. Fig. 3 (right), agrees well with a similar calcu- 
lation performed in the NIL model [15], as was expected since the latter includes quark degrees of 
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freedom in the vacuum by construction. In this case, the mentioned new phase - where the quarks 
are aheady deconfined, but chiral symmetry is still broken - should also occur. Recent calculations 
in a modified version of the NJL model suggest the existence of a quantitatively new, third, scenario 
where the gap between the two transitions is very small, while the nearly-common critical temper- 
ature increases with the magnetic field [16]. Lattice simulations of QCD with dynamical quarks 
also indicate a rise of the critical temperatures and show no evidence of the splitting [17]. However, 
in this lattice simulation the dynamical quarks are still heavier compared to realistic values. On the 
other hand, a holographic calculation points out to the existence of a small splitting between the 
transitions [14]. 
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